     


figure(10);
 count = 0; 
 for h=[1 4 8 12]
    count = count+1;
 subplot(2,2,count);

plot([ones(size(MSFE_var{h}+h,2),1) prctile(MSFE_var{h}./MSFE_ar{h},[16 50 84],1)'],'r');  datetick('x','keepticks')
anni=datevec(nn(16+1:80:end,:));
set(gca, 'XTick', 16+1:80:length(nn));
set(gca, 'XTickLabel', anni(:,1));  
 
 hold on
 plot([ones(size(MSFE_var{h}+h,2),1) prctile(MSFE_var{h}./MSFE_ar{h},[50 ],1)'],'k'); datetick('x','keepticks') 
anni=datevec(nn(16+1:80:end,:));
set(gca, 'XTick', 16+1:80:length(nn));
set(gca, 'XTickLabel', anni(:,1));  

 % hold on
% plot([datenum(1914,9,15) datenum(1914,9,15)], [0 9]); axis([696410 732021 0 9] )
% plot([datenum(1931,6,15) datenum(1931,6,15)], [0 9]); axis([696410 732021 0 9] )
% plot([datenum(1941,9,15) datenum(1941,9,15)], [0 9]); axis([696410 732021 0 9] )
% plot([datenum(1951,3,15) datenum(1951,3,15)], [0 9]); axis([696410 732021 0 9] )
% plot([datenum(1971,9,15) datenum(1971,9,15)], [0 9]); axis([696410 732021 0 9] )
% plot([datenum(1983,12,15) datenum(1983,12,15)], [0 9]); axis([696410 732021 0 9] )
% hold off
% title('12-QUARTERS HORIZON')
% % set(gcf,'DefaultAxesColorOrder',[0 0 0;1 0 0;1 0 0;0 0 1]); 
end
 